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ABSTRACT 



Context. With the new generation of high-resolution X-ray spectrometers the understanding of warm absorbers in active galactic 
nuclei has improved considerably. However, the important questions regarding the distance and structure of the photoionised wind 
remain unsolved. 

Aims. To constrain the distance of the photoionised wind, we study the variability of the continuum, absorption, and, emission in one 
of the brightest and most variable low-luminosity AGN: the narrow-line Seyfert galaxy NGC 4051. 

Methods. We analyse two observations taken with the Low Energy Transmission Grating Spectrometer of Chandra. We investigated 
the spectral response to a sudden flux decrease by a factor of 5, which occurred during the second observation. 
Results. We detect a highly ionised absorption component with an outflow velocity of -4670 kms"', one of the highest outflow 
velocity components observed in a Seyfert 1 galaxy. Furthermore, this is one of the only observations whereby the X-ray observed 
absorption component is unaccompanied by a corresponding UV absorption component with the same outflow velocity. The spectra 
contain a relativistic O vm Lyo- line with properties similar to those determined for this source with XMM-Newton, and four absorption 
components spanning a range in ionisation parameter f between 0.07 and 3.19 (log values, and units of 10"' W m). An emission 
component producing radiative recombination continua of C vi and C v appears during the low state. The black body temperature 
decreases with the drop in flux observed in the second observation. 

Conclusions. For all absorber components we exclude that the ionisation parameter linearly responded to the decrease in flux by 
a factor of 5. The variability of the absorber suggest that at least three out of four detected components are located in the range 
0.02 - 1 pc. For one component we only have a lower limit of 0.3 pc. These distances are different from earlier suggestions. 

Key words, galaxies: Seyfert - quasars; individual; NGC 4051 - galaxies; active - X-rays; galaxies 



1. Introduction 

The stupendous amount of energy emitted by an active galactic 
nucleus (AGN) is released by gas that flows towards a super- 
massive black hole in the centre of a galaxy, presumably via 
an accretion disc. The average energy of the emitted radiation 
increases with decreasing distance to the super-massive black 
hole, thus X-rays provide the best probe of the gas flow in the 
immediate surroundings of the black hole. From the presence of 
jets, as well as from blue-shifted ultraviolet absorption lines, we 
learn that there is gas flowing away from the black hole, in ad- 
dition to the flow towards it. With the availability of X-ray spec- 
trographs onboard Chandra and XMM-Newton, we can study 
the lines emitted in the X-ray band from these outflows close to 
the black hole. The presence of multiple absorption line systems 
that differ in their level of ionisation or in their outflow velocity 
constrains th e geometry of the outflo w (for exam ple NGC 3783, 
iNetzer et al.ll2003; NGC 5548. Steenbrugge et al. 2005) . 

In this paper we describe the X-ray spectra of the narrow- 
line Seyfert 1 galaxy NGC 4051, obtained with the Low-Energy 
Transmission Grating (LETG) on board of the Chandra satel- 
lite. Em i ssion from the nucleus of NGC 405 1 was detected by 
iHubbld (Il926h . and the galaxy is o ne of the 12 discussed in 
the seminal paper by ISevfera ( Il943h . As a member of the Ursa 



Major cluster (jT ull y et al.lll996l) . it has a distance of 18.6 Mpc 
(Tullv & Pierce 2000) and a heliocentric velocity of 700 kms"' 
( Verheiien & Sancisi.2001.) . NGC 4051 is a relatively bright op- 
tical object, with V ^ 13.5. The soft excess observed with the 
Reflection Grating spectrometers (RGS) onboard XMNi-Newton 
is not described well by a multi-temperature disc: Ogle et aT] 
(|2004 concludes that the soft excess is due to emission from 
relativistically broadened O vm emission lines and its related ra- 
diative recombination continuum (RRC). However. lPounds et al.l 
('200?) fitted the EPIC data of the same observation satisfactorily 
with a power law plus two black body components. 

Two absorption line systems, at -2340 + 130 and -600 + 
130 kms were detected in the X-ray spec trum with the 
High Energy Transmission Grating of Chandra dCoIIinge et al.l 
2001). Nine absorption systems with a small range in out- 
flow velocity are detected in the ultraviolet with the Space 
Telescope Imaging Spectrograph (STIS) and the Far Ultraviolet 
Spectroscopic Explorer (FUSE). One of the absorbers observed 
in the UV may coincide with the -600 + 130 kms"' X-ray ab- 
sorption system, all the others have s maller outflow velocities 
(ICoIIinge et aI.ll200UlKaspietani2004i) . 

In a prelimina r y ana lysis of our first LETGS spectrum. 



nar y 

Ivan der Meer et alj (l2003h did not detect the -2340 km s ' ve- 
locity component, but found indications for the presence of an 
outflow component at an even higher velocity of -4500 km s"' . 
The second observation, taken when NGC 405 1 was in a much 
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Fig. 1. The light curve in the zeroth spectral order (top panels) and hardness ratio (bottom panels) for the observations of Jan 2002 
and July 2003. The hardness ratio is defined as the ratio of the 2-8 and 8 - 32 A counts. The top panel has a bin size of 350 s and 
the bottom panel has a bin size of 2000 s. Dotted lines denote the different flux states and define the intervals used in creating the 
four spectra. 



higher flux state, allows us to investigate the ionised outflow both 
in terms of outflow velocity and ionisation structure. In this pa- 
per we report the analysis of both observations with the LETGS. 

2. Observations and data reduction 

In this paper we describe the t wo observations of NGC 4051 
taken with the LETGS (Bri nkman et all |2000|) onboard of 
Chandra. The detector used was the High-Resolution Camera 
(HRC-S). The observations span the time periods December 
31, 2001 17:40:41 to January 1, 2002 19:51:46 and July 23, 
2003 00:34:15 to July 24 03:24:43, and have a net exposure 
time of 94.2 and 96.6 ks, respectively. Henceforth we will we 
refer to these observations as the Jan 2002 and Jul 2003 ob- 
servations. The data reduction method is described in detail by 
iKaastra et alj ( l2002h . After event selection we produced back- 
ground subtracted spectra and light curves. As the HRC-S de- 
tector has limited energy resolution, higher spectral orders can- 
not be separated in the data but are taken into account during 
spectral fitting. We only use spectra between 2 A and 80 A, as 
the Galactic absorption attenuates the spectrum strongly above 
~ 80 A. 

3. Variability 

Fig. [T] shows the light curves for each observation. For the Jan 
2002 observation, the average brightness of the source increased 
noticeably in the second half of the exposure. The source was 
brighter still during the first ~ 65 ks of the Jul 2003 observation 
and then dropped in about 3 ks to a lower flux level, which is 
similar to the flux level of the the second half of the first obser- 
vation. The flux recovered gradually to the high level of the first 
half of the observation during the final 20 ks. On top of these 
long-term variations we also detect variability with smaller, but 
significant, amplitude on timescales of a few thousand seconds. 
These rapid variations have smaller amplitudes after the rapid 
flux decrease of the Jul 2003 observation. In the analysis below 
we will discuss the spectra at different flux levels, defined by the 
lightcurve of Fig. [1] A and B are the first and second half of the 
Jan 2002 observation; C and D are the first and second part of 
the Jul 2003 observation. 



In Fig.[T]we also show the hardness ratio, defined as the ratio 
between the count rates in the 2-8 A and 8-32 A band passes. The 
hardness variation between different flux levels is small, with 
HR^O.239 ± 0.14, 0.203 + 0.008, 0.161 ± 0.003, 0.134 ± 0.021 
for intervals A, B, C and D, respectively with their Icr error. 

4. Spectral data analysis and results 

From the above (Fig.[TJ it is clear that NGC 4051 shows signifi- 
cant luminosity variability during our observations. As the level 
of the variations is quite strong (e.g., a factor of five in flux be- 
tween parts C and D of the July 2003 observation), it is not useful 
to analyse the integrated spectra as a whole, but we will rather 
analyse spectra A, B, C and D separately, with a focus on the 
spectral differences between these parts. Statistics prohibit us to 
do spectral analysis down to the shortest variability timescales 
of less than 1000 s. 

For all spectral fitting we used the SPEX code developed and 
regularly updated with new models and atomic data at SRON 
dKaastra et aLi[l996 !Pl. As the flux is low in some of the spec- 
tra, we cannot use minimisation, instead we minimise the C- 
statistic. Note that that the uncertainty in absolute calibration for 
the Chandra LETGS is 13 mA. Therefore a velocity difference 
smaller than ~ 200 km s"' between spectra A/B and C/D is un- 
likely due to an intrinsic velocity shift. 

For our spectral model we included different components: a 
broad-band continuum, narrow and broad emission features, and 
photoionised absorption. The need for these components and the 
resulting parameters will be discussed separately. 

Noting that spectrum C has the highest signal-to-noise of the 
four spectra, we decided to fit this spectrum first and then use the 
best fit model as a template for fitting spectrum B. The signal-to- 
noise ratios for spectra A and D are poor, therefore we only left 
a limited number of parameters free when fitting these spectra. 

The best fit parameters for our spectral model are listed in 
Table[T] All models were corrected for Gal actic absorption with 
acolumndensity of 1.31x10^"* mT^ jElvis e t al. 1989) and for the 
cosmological redshift (including peculiar motion) of 700 km s"' . 
All errors correspond to AC = 1 (68 % confidence). The best fit 

' http://www.sron.nl/spex 
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Table 1. The best fit parameters for the different spectral components for the four spectra. All wavelengths and velocities refer to 
the rest frame of NGC 4051. Numbers in square brackets are calculated from the fit parameters; numbers in round brackets were 
kept fixed during the fit. 





A 


B 


C 


D 


Continuum 










• power law 










Flux (2-10 keV)° 


0.83 + 0.07 


1.36 ±0.08 


2.93 ± 0.09 


0.48 ± 0.05 


Photon index T 


1.65 + 0.07 


2.00 ± 0.05 


2.24 ± 0.02 


2.44 ± 0.05 


• modified black body 










Flux (0.2 - 10 keVf 


0.49 ± 0.08 


0.66 ±0.14 


2.0 ± 0.2 


0.37 ± 0.09 


kT (eV) 


95 ±6 


96 ±8 


143 ±4 


90 ±5 


Absorption components 










• warm absorber 1 










Nu (10-^ 


- 


- 


0.12 ±0.04 


(0.12) 


log^ (10-' Wm) 


- 


- 


0.07 ±0.13 


0.3 ±0.2 


Vturb (kms"') 


- 


- 


300 ± 60 


(300) 


Vout (kms-') 


- 


- 


-210 ±70 


(-210) 


• warm absorber 2 










Nu (10^5 m-2) 


(0.17) 


0.17 ±0.05 


0.29 ± 0.05 


(0.29) 


log^ (10-' Wm) 


0.95 + 0.25 


0.49 ± 0.20 


0.87 ± 0.09 


0.52 ±0.13 


Vturb (km s-') 


(230) 


230 ± 120 


110 ±20 


(110) 


Vout (kms"') 


(-330) 


-330 ± 140 


-200 ± 30 


(-200) 


• warm absorber 3 










Nu (10^5 m-2) 


(0.6) 


0.6 ±0.3 


0.8 ±0.4 


(0.8) 


logf (10-' Wm) 


2.23 + 0.17 


1.98 ±0.22 


2.32 ±0.13 


7 4.+0.6 


Vturb (kms-') 


(100) 


100 ± 50 


90 ±30 


(90) 


Vout (kms"') 


(-610) 


-610± 110 


-580 ± 50 


(-580) 


• warm absorber 4 










Nu (W^ m ^) 


(9) 


9 ± 5 


20 ± 10 


(20) 


log ^(10"' Wm) 


3.06 ±0.10 


2.92 ± 0.17 


3.19 ±0.09 


3.2 ± 0.2 


Vturb (kms ') 


(120) 


120 ± 80 


19 ± 13 


(19) 


Vout (kms ') 


(-4590) 


-4590 ±180 


-4670 ± 150 


(-4670) 


Emission lines 










• (J VIII Lya { relativistic) , 


c 








rlux 


iZ + Z 


09 -1- ^ 

ZZ ± J 




i.j ± 1.1 


cj (index Laor profile) 




_i_ n ^ 

J.Z ± U.J 


Z.UD ± U.ZO 




FWHM (A) 


[10] 


[10] 


[3.0] 


[3.0] 


* 0\m Lya (BLR f 










Flux"* 


0.1 + 0.2 


U.O ± O.Z 


i.8 ± U.6 


U.S ± U.4 


FWHM (A) 


(0.45) 


(0.45) 


0.45 ± 0.14 


(0.45) 


• Ovii triplet (BLR) 










Flux'' 


1.4 + 0.6 


0.0 ± 0.7 


6.2 ± 1.4 


2.2 ± 1.2 


FWHM (A) 


(1.7) 


(1.7) 


1.7 ±0.4 


(1.7) 


^(A) 


(22.11) 


(22.11) 


22.11 ±0.13 


(22.11) 


• unforbidden (narrow) 










Flux'' 


0.8 ± 0.2 


0.9 ± 0.2 


1.0 ±0.2 


1.3 ±0.3 


A (A) 22.087 ± 0.005 


22.081 ±0.006 


22.090 ± 0.009 


(22.090) 


• Cm Lya (BLR)' 










Flux'' 


0.6 ± 0.4 


1.0 ±0.5 


2.4 ±0.8 


3.1 ±0.7 


FWHM (A) 


(0.36) 


(0.36) 


0.36 ±0.13 


(0.36) 



" Absorption-corrected flux in 10 Wm ^. 
* Wavelength is frozen to the rest wavelength of 18.967 A (for O vm Lycf). 
^ And an inclination angle of 48°. 
'' Absorption-corrected photon flux in photons m"^ s"'. 
'' Wavelength is frozen to the rest wavelength of 33.736 A. 



values for the C-statistics are 2700 (2476), 2658 (2466), 2738 
(2599), and 2870 (2610) for components A-D; the numbers in 
brackets are the degrees of freedom. In Table|2]we show the in- 
crease AC of the C-statistic when a given spectral component is 
omitted from the model (and re-fitting the spectrum; the omis- 
sion of a given component can sometimes be partly compensated 
by adjusting the parameters of the other components). We did 
this for spectra B and C, the best quality spectra. 



Finally, in spectrum D we afterwards added a spectral com- 
ponent fitting the radiative recombination continua (RRCs). 
Only the C v and C vi RRCs are present, and can be identified 
by eye in Fig.|2] The temperature and emission measures n^nxV 
with V the emitting volume, the electron density, nx the den- 
sity of the recombining ion X giving rise to the observed transi- 
tions are given in Table [3] 
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Fig. 2. Flux (in units of photons nT^ s"' A"') and our best spectral model. All spectra are binned by a factor of 2 for clarity of 
display. Note the different continuum levels and shapes between the spectra, especially around the upturn at ~ 16 A. In spectrum D 
C VI and weaker C v RRCs at 25.43 A and 31.80 A respectively, are observed. 
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Table 2. Significances of the different spectral components. We 
list the increase in the C-statistic when the component is omitted 
from the model and the spectrum is refitted. 



Component 


Spectrum B 


Spectrum C 


Modified black body 


7.36 


88.97 


rel. VIII 


16.84 


11.32 


BLR vm 





16.69 


BLR VII 





31.49 


narrow vii forbidden 


36.08 


24.60 


BLR C VI 


4.07 


12.01 


warm absorber 1 




14.53 


warm absorber 2 


12.74 


17.83 


warm absorber 3 


12.75 


14.55 


warm absorber 4 


45.19 


32.89 



Table 3. The best fit parameters for the RRCs in spectrum D. 
The velocity refers to the rest frame of NGC 405 1 . 



Parameter 

kT (eV) 

CV 

Cvi" 

Nvi" 

Nvn" 

Ovn" 

Ovm" 

velocity (km s"') 



value 
5±2 
210 ± 110 
390 ± 80 
40 ±60 
30 ±30 
0±20 
10± 14 
+ 1300 ±200 



" Emission measures in 10*'* m ^, e.g., X=C vii produces the RRC 
observed for C vi. 



4.1. Continuum 

In spectrum C the continuum cannot be fitted with only a 
power law. The soft excess can be modelled by adding a black 
body modified by coherent Compton scattering (iKaastra & Bard 
11989 ). Examining all spectra, we find that the power-law pho- 
ton index becomes softer as time progresses, and is thus not a 
simple function of luminosity. The modified black body temper- 
ature seems to track the luminosity variations; the temperature is 
highest for the high state and drops back to a lower value after 
the drop in flux just before the start of spectrum D. 



4.2. The ionised absorber 

Thanks to the excellent spectral resolution of the LETG instru- 
ment we can study the warm absorber in detail. We fitted the 
warm absorber with the xabs model in SPEX. In this model, the 
ionic column densities are not independent quantities, but are 
linked via the ionisation parameter ^ = L/nr^, where L is the 
source luminosity, n the hydrogen density and r is the distance 
from the ionising source. For the spectral energy distribution of 
the ionising source we used the standard one included in SPEX 
(see dot-dash line of figure 1 in Steenbrugge et al. 2005). In ad- 
dition to the ionisation parameter, the fit parameters are the hy- 
drogen column density of the absorber Nh, the outflow velocity 
Voui, and the turbulent velocity Vmrb- The advantage of the xabs 
model is that all relevant ions are taken into account including 
also the ones with the weakest absorption features. In the fit the 
abundances were left froz en to the solar abundances given by 
lAnders & Grevessd(ll989l) . 




22 23 
Wavelength (A) 



Fig. 3. Detail of spectrum C fitted with 3 warm absorber compo- 
nents (thin line) and 4 warm absorber components (thick line). 
Binning as in Fig.|2] The best fit 3 component model does under- 
predict the depth of the O v and O iv lines. 



Table [T] lists the best fit values for the four different spec- 
tra. To adequately describe the warm absorber in spectrum C we 
need 4 components. The ionisation parameters range between 
0.07 and 3.19 (logarithmic values in 10"'' W m, used through- 
out this paper). Adding a 5''' absorption component for spectrum 
C does not significantly improve the fit, namely A C decreases 
by only 1.91. Fig. [3] compares the best fit model with 3 and 4 
absorption components and Table |4]gives the best fit parameters 
for the three-absorber-component fit to spectrum C. The fit with 
4 warm absorber components better fits the O v and O iv absorp- 
tion lines. In the remaining three spectra due to the lower signal 
to noise ratios, not all four absorption components are detected. 



Table 4. The best fit parameters for a three-absorber-component 
model for spectrum C. Components 1-3 correspond to compo- 
nents 2 - 4 in Tables [1] and |2] 



Parameter value 

• power law 

Flux (2 - 10 keV)° 2.86 ± 0.04 

Photon index T 2.21 ± 0.02 

• modified black body 

flux (0.2 - 10 keV)° 1.31 ± 0.14 
kT (eV) 140 ± 4 



Absorption components 

• warm absorber 1 



Nh (10^5 m-2) 
log^ (10"' Wm) 
VtLirb (kms-') 
Vout (kms-') 

• warm absorber 2 
Nh (10^^ m-2) 
log^ (10-' Wm) 
vtmb (kms"') 

Vo„t (kms"') 

• warm absorber 3 
Nh (10^^ m-2) 
log^(10-' Wm) 
Vturb (kms"') 

Vout (kms-') 



0.35 ±0.14 
0.73 ± 0.06 
147 ± 21 
-213 ± 27 



0.72 : 
2.27 ± 

100: 

-548 

23 ± 
3.13 ± 



tO.3 

:0.06 

t 28 
±45 

7.5 

:0.04 



-4530 ± 80 



' Absorption-corrected flux in 10 Wm 
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For spectrum B, we need three components. By comparing 
the parameters of these components to the parameters of the 
components for spectrum C, we find that component 1, which 
has the smallest hydrogen column density, is lacking in spec- 
trum B. 

The statistical quality of spectra A and D are much poorer; 
in fitting these spectra, we fixed several parameters to the best-fit 
values obtained for spectra B and C, respectively. It is quite rea- 
sonable to keep the outflow and turbulent velocity constant be- 
tween observations A/B and also between C/D. When we leave 
both the column density and the ionisation parameter free in 
spectrum A, only component 2 is marginally detected. These 
non-detections do not imply that the warm absorber has disap- 
peared, but merely that the statistical quality of spectra A and 
D are insufficient to determine its parameters with accuracy. For 
this reason, we make the additional assumption that the column 
density of the absorber does not change. Any changes however, 
are likely to be in the ionization parameter (that is, if it can re- 
spond fast enough to the continuum flux variations). 

Comparing spectra C and D, only for component 2 a 2.3cr 
drop in ionisation parameter is found; the other components are 
consistent with no change in the warm absorber Comparing 
spectra A and B, we find only marginally higher ionisation pa- 
rameters for components 2 - 4 in spectrum A. 

Fig. ID shows the change in C-statistics for spectra A and D 
when the ionisation parameter for each component is separately 
varied. Step size is 0.1 in log ^ and at each step the goodness 
of fit is calculated. It is clear from this figure that no significant 
changes can be detected for components 1, 3 and 4. We discuss 
the variability of component 2 in more detail later. 

In Fig.|5]we present the model absorption spectrum for each 
component, labelling the deepest lines. It thus allows the reader 
to see which component contributes to which line. The contin- 
uum was renormalised to ease identification. 

4.3. Broad emission iines 

Besides absorption lines, strong emission features are present in 
the four spectra. We chose two different methods for fitting this 
excess emission: a model with broadened emission lines and a 
model in which the absorbers only partially cover the contin- 
uum. We first describe the model with broadened emission lines. 
As was the case for fitting the absorber, we first fitted spectrum C 
and then used its best fit values as a starting point to fit the other 
spectra (see Table [T]for best fit parameter values). The presence 
of a broad O vm Lya emission feature around 19 A is illustrated 
in Fig. |6] which shows a detail of the spectra B and C includ- 
ing the best fit model with and without relativistically broadened 
O VIII Lyff. Fig. |2] shows the model emission spectrum for spec- 
tra B and C. Excess emission with at least a 1 A width, due to 
the O VII triplet near 22 A is present in spectra A, C and D. In 
addition, there is significant excess emission in spectra C and D 
near 34 A, the region of the ls-2p transition in C vi. 

We modelled the excess emission at 1 9 A (see F ig.!?} with a 
relativistically broadened hne, fol lowing lOgle et al.l t2004) and 
iBranduardi-Ravmont et al.l (12001 1) . The model used is a narrow 
emission lin e convolved with the relativistic disc line profile of 
lLao3 ([1991 . The inner (ri) and outer (r2) disc radii were fixed 
in our models to the default values (ri - 1 .234GM / and r2 - 
400GM/c2). The rest-frame wavelength of the line was frozen to 
18.969 A, the wavelength of O vm Lya. 

We have tried fits with different disc inclinations /, but in 
all fits the inclination converged to 48°, within a few degrees. 




component 1 
• component 2 
A component 3 
■ component 4 




log 5 

Fig. 4. Change (AC) in C-statistic for spectrum A (top panel) 
and spectrum D (bottom panel) when the ionisation parameter ^ 
of each component is varied. Open circles indicate the reference 
values of each component (for spectrum D, the value of ^ is the 
one determined from spectrum C; for spectrum A, the value of ^ 
is the one determined from spectrum B). 



This value is consistent with the one found byl Ogleet al.](l2004 . 
The disc inclination is not expected to change significantly over 
a couple of years. For that reason, and in order to reduce the 
number of free parameters, we have kept the inclination frozen 
to that value in all presented fits. The disc emissivity slope q 
(with R the distance to the black hole) was allowed to vary, 
however. 

In addition to a relativistic Ovm Lya line, we added a 
Gaussian-shaped broad O vm line Lya line (FWHM 0.45 A) to 
our model for spectrum C and D. Further broad lines correspond- 
ing to the O vii triplet and C vi Lya line were also needed. For 
C VI Lya, we have frozen the wavelength to the rest wavelength 
of 33.736 A. For Ovii we kept the wavelength as a free pa- 
rameter: the line is a blend of the O vii triplet, and the relative 
contributions of these lines depend on the physical conditions 
in the source. Such broadened emission lines were first detected 
in NGC 5548 by Kaastra et a l. (2() 02), and confirmed for that 
source by Steenbrugge et alj (l2005h . and were detected in this 
source bv .Ogle et al.i (i2004ft . 



K. C. Steenbrugge et al.: High-resolution X-ray spectroscopy of NGC 4051 



in 



Z3 in 

o 

O 



o 

CO 



o 

C\J 



8 ° 




15 



20 

Wavelength (A) 



25 



Fig. 6. Detail of spectra B and C with the best fit model (thick line) and the best fit model without the O viii Lya relativistically 
broadened emission line (thin line). 











absorber 1 1 
absorber 2 






r n 






absorber 3 


M 






















pi 


T 


1 


1 


IP 




' r 


r ' 








■| 


■ 


absor 


1 

e 


|| 

4 

1 


|1 


n 


1 




I 

O 


1 




1 


j 

1 





10 20 30 

Restframe wavelength (A) 



Fig. 5. The model absorption spectra for each of the 4 different 
absorption components as detected in spectrum C. The contin- 
uum was modelled to be flat and was renormalised per compo- 
nent. 



4.4. Partial covering absorption model 

As an alternative to fitting extra emission lines, we fitted the 
spectra assuming that the absorber does not fully cover the con- 
tinuum emission. Considering that the excess emission is most 
significantly detected in spectrum B and C, we first fitted these 
two spectra (see Table [T] and Fig. |6]l. In this fit there are no 
broad emission lines fitted. We fixed the covering factor for the 
3 or 4 detected absorption components (spectrum B and C, re- 
spectively) to 0.75, and then refitted the continuum parameters 
(power law and modified black body) and the hydrogen column 
density and ionization parameter of the detected absorber com- 
ponents. Afterwards we also allowed the covering factor for each 
absorption component free to vary in the fit. 

For spectrum C the best fit covering factors for three of the 
four absorption components is 1 . For warm absorber component 
2 the best fit covering factor is 0.7 + 0.13, with an increased 
hydrogen column density of 0.63 x 10^^ m^- (compared with 
0.29 X lO^*" m"^ for a covering factor of 1). Forcing all cover- 
ing factors to 1 increases AC by 17. However, the best fit with 
partial covering increases AC by 70 for a decrease of 7 in de- 
grees of freedom, compared with the best fit model with broad 
emission lines (4 extra free parameters in the partial absorption 
fit compared to 1 1 in the model with broad emission lines). A 
comparison between the best fit with broad emission lines and a 
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Fig. 8. Top: detail of the spectrum for the best fit (thick line) model with broad emission lines and the best fit (thin lines) assuming 
the absorber has a covering factor of .75 for spectrum B. Bottom: same as above but for spectrum C. 



fit with all absorbers having a covering factor of 0.75 can be seen 
in the bottom panel of Fig. [8] Note that the continuum long-ward 
of ~24 A is consistently over-estimated in the fit assuming that 
the absorber only covers part of the continuum emission. 

For spectrum B we followed the same procedure, but here 
the best fit for the covering factors of the three absorption com- 
ponents detected are 1 , even if we redid the fitting starting from a 
covering factor of only 0.25. The best fit without broad emission 
lines increases AC by 21 and is shown in Fig.|6l Forcing the three 
absorbers to have a partial covering factor of 0.75 worsens the fit 
by a further 35 in AC. This last fit is compared with the best fit 
with broad emission lines in the upper panel of Fig. [8] Note that 
again the continuum is over-estimated at long wavelengths, fur- 
thermore the model assuming partial covering cannot reproduce 
the steep rise at 16 A. 

4.5. Narrow emission lines 

The Ovii forbidden line at 22.10 A is detected with at least 
4cr significance in each individual spectrum. There is no indi- 
cation for significant flux variability or velocity broadening for 
this line. 

In spectrum D a strong, ~ 5cr significant, RRC from C vi 
(rest wavelength 25.30 A) and a ~ 2cr significant RRC of Cv 
(rest wavelength 31.62 A) are present (see Fig. |2] and Table |3]l. 



In addition, we have determined limits to the strength of other 
expected RRCs of oxygen and nitrogen ions (see Table O. The 
temperature of the recombining photoionised gas is low but not 
well determined (T - 5 + 2 eV). In our final fit we fix the tem- 
perature to 5 eV, in order to better assess the significance of this 
component. There is a strong correlation between the adopted 
temperature and the strength of the emission measure. The RRCs 
show a redshift with respect to the rest frame of the galaxy of 
~ 1300 kms'. 

5. Discussion 

5.1. Continuum 

The uniqueness of the transition between part C and D is that the 
rather abrupt flux change is preceded and followed by long peri- 
ods of relative quiescence. This allows us to study the response 
of the continuum, the broad and narrow emission lines, and the 
warm absorber to this sudden flux drop. We first describe the dif- 
ference in the continuum. Fig. |9] shows the difference in overall 
continuum profile for spectrum C and D. 

Between spectrum C and D the temperature measured for 
the modified black body significantly decreases, and the power- 
law slope becomes softer The temperature remained constant 
and the power-law slope softened between spectrum A and B. 
For Seyfert 1 galaxies, normally, the spectral slope hardens 
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Fig. 7. Model emission spectra for spectrum B and C, without 
Galactic absorption and the warm absorber components. Dashed 
line: power law continuum; dotted line: modified black body 
continuum; dash-dotted line: sum of power law and modified 
black body continuum; thin solid lines: relativistic, broad and 
narrow emission line contributions; thick solid line: total unab- 
sorbed emission spectrum. 



with decreasing l uminosity with little or no change in the black 
body component. [Pounds et al.l (l2004t) studying two EPIC spec- 
tra taken in a high and a low state of NGC 405 1 found that the 
higher temperature black body component disappeared in the 
low state. The authors found that the power-law slope did not 
vary between the high and low state, which is different from our 
results. 

Changes in the warm absorber properties could mimic a tem- 
perature change of the modified black body component. To en- 
sure that in the spectra there really is a temperature change we 
did the following tests. We refitted spectrum D, but freezing the 
temperature of the modified black body to the best fit value of 
spectrum C. In the fit the column density and ionisation parame- 
ter of the four absorption components were left free to vary. The 
fit thus obtained was substantially worse (increasing C-statistics 
by 102), and the fitted ionisation parameters were very similar 
to those obtained in the original fit of spectrum D. Afterwards, 
we lowered the ionisation parameter of the only significantly de- 
tected warm absorber component in spectrum D by a (linear) 
factor 5. This decrease matches the decrease in luminosity ob- 
served. Leaving the ionisation parameters of the absorber and 




Fig. 9. The spectrum for part C and D, clearly showing the diff'er- 
ence in continuum shape. The difference is due to a lower flux in 
the power-law component and a lower temperature for the mod- 
ified black body component in spectrum D. 



the temperature of the modified black body free in the fit,we 
ag ain derive our origi nal fit parameters given in Table [1] 

lOgle et al.1 (12004') studying the XMM-Newton RGS spec- 
trum, find that a multi-temperature disc model is a poor fit to 
the soft excess, and that instead it is well fitted by the relativis- 
tically broadened O viii line series. We tested this model on our 
highest signal-to-noise spectrum, C, by adding to the O vin Lya 
line also the Ly/3, Lyy and the RRC emission lines relativisti- 
cally broadened. As none of these emission features are obvious 
from a visual inspection, we fixed the relative normalisation of 
the different lines to the Lya component: 0.25 for Ly/3, 0. 1 for 
Lyy and 0.21 for the RRC. The parameters q and / were free in 
the fit, but the same for the four different emission lines. The fit 
obtained was significantly worse (increasing C-statistics by 195) 
for the same number of parameters than a fit with a modified 
black body. 

We conclude that there is emission from a modified black 
body present in the LETGS spectrum, and that the temperature 
of this component drops with decreasing luminosity between 
spectrum C and D. 



5.2. Relativistically broadened Oviii Lya line 

The relativistically broadene d Oviii Lyg line which explains 
part of the RGS soft excess dOgle et al.l 12004) is significantly 
detected in all spectra (see Table [1). The measured line fluxes 
do not follow the power law or modified black body continuum 
flux changes. The emissivity q, assuming R"^ , is similar to the 
value q2 -5.0 measured bv lOgle et al.l (|2004|) for spectrum B, 
but not for spectrum C. We note that for spectrum A and D the 
emissi vity was frozen to the value for spectrum B and C respec- 
tively. 'O gle et all (|2004|) in detail discusses the steep emissivity 
derived as well as the physical parameters in the accretion disc 
derived from the line profile and the lack of C vi Lya relativistic 
line. The latter places a lower Umit on the ionisation parameter 
Summarized, the relativistic Ovm Lya is emitted between L6 
and 6 gravitational radii, ^ > 3.4 and ng < lO^'' m"-'. 

The inclination angle for all spectra was cons istent with 
48°; very similar to the inclination measured by lOgle et all 
(i2004i) . The inclination angle for the nan^ow line region (NLR) 
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as derived by IChristopoulou et aTl (Il997h from narrow band 
[O m] imaging and long-slit spectroscopy using the Manchester 
Echelle Spectrograph is about 50°. This indicates that the nar- 
row line region, and the gas producing the relativistically broad- 
ened emission line lie along the same line of sight. The narrow 
line region is therefore possibly an extension of the gas produc- 
ing the relativistically broadened emission lines. If this is the 
case, then the gas forming the relativistically broadened emis- 
sion line is outflowin g, similar to the absorptio n components. 
From the kinematics, Christopoulou et aTl (1 19971) conclude that 
the imaged [Oiii] gas is outflowin g. The inclination angle o f 
the host galaxy is at medium angles ( MoUenhofi' & H eidt 2001). 
thus gas observed from the relativistic and narrow emission lines 
are closely aligned with the host galaxy. 

Relati vistic line profiles in the soft X-r ay band were first re- 
ported bv iBranduardi-Ravmont et al] (1200 1.) in the narrow lines 
Seyfert 1 galaxies MCG -6-30-15 and Mrk 766. Interestingly, 
the relativistically broadened Oviii, Nvii and Cvi Lya lines 
have only been convincingly detected in narrow line Seyfert 1 
galaxies. For MCG -6-30-15 and Mrk 766 the disc line param- 
eters derived from the soft X-ray band and the Fe-K line are in 
good agreement. The XMM-Newton obse rvation of NGC 4 05 1 
that showed a strong relativistic O viii line (Ogle et al. 2004) did 
not show a significant relativistic Fe-K line (Pounds et al. 2004). 
Unfortunately the LETGS is not sensitive enough in the the Fe-K 
band to directly compare the disc geometry as derived from iron 
and oxygen lines. For the physical implications of th e relativis- 
tic line s in NGC 405 1 we refer to the discussion by lOgle et aTl 
(12004 . 



5.3. Broad emission lines 

The broad O vii line, centred on the forbidden line, is detected 
in spect rum C, bu t the l ine is weaker or absent in the other 
spectra. ,Ogle et al.l (|2004|) detected a broadened O vii emission 
line with a fuU width half maximum (FWHM) of 1 1000 + 3000 
km s"'. This is consistent with the FWHM we derive: 24000 
+ 6000 km s"'. Both values are rather larger than the FWHM 
dete rmined from the HyS and He ii broad lines measured in the 
UV ("Pet erson et ani2000 ^. If the broad line region (BLR) has a 
stratified ionisation structure, then it is expected that the higher 
ionised O vii line will have a broader width than either the H/3 
or lowly ionised He ii. At 1 .7 A the width of the broadened O vii 
line exceeds the separation (0.5 A) between the forbidden and 
resonance line. The broadened resonance (and intercombination) 
lines contribute to the flux and hence the width measured. 

We also detect a broadened C vi line with a width of 3200 + 
1000 km s which is smaller by a factor of ~7 than the Ovii 
line. This value is fairly similar to the FWHM derived for the 
H/3 line (1100 + 190, Peterson et al. 2000) and the wi dth mea- 
sure d for this line in the RGS spectrum (1200 km s '. IOgle et alj 
l2004i) . The smaller width is consistent in terms of a stratified 
BLR model, with the fact that C vi is formed at lower ionisation 
parameters. 

In spectrum C on top of the relativistically broadened O vm 
Lya emission line there is a broadened O viii Lye emission line. 

ICostantini eTal] (|2007|) sfiidying Mrk 279 did detect broad- 
ened emission lines with better statistics than is the case for 
NGC 4051. The authors conclude that the broad X-ray emis- 
sion lines originate from a stratified BLR. The above findings 
for NGC 405 1 are consistent with this conclusion. 



5.4. Narrow emission features 

During our observations we detect various narrow emission fea- 
tures in the spectra. The narrow forbidden O vn line was signif- 
icantly detected in all spectra, and all flux values are consistent 
with each other to within Icr (see TablelTJ. This is consistent with 
the picture that the O vii forbidden line is formed at the narrow 
line region further out. The width of this line (it is unresolved), 
precludes it being formed as part of the outflow. 

In spectrum D we detect the C vi RRC and the C v RRC. In 
the other three spectra these features are not detected. The shapes 
and centroids of individual RRCs are determined by the tempera- 
ture of the recombining gas, the emission measure of the relevant 
ion, and by the motion (velocity broadening) of the recombining 
gas. We derived the temperature and emission measure for the 
RRCs assuming no velocity broadening. We have investigated 
the eff'ect of broadening by including a velocity broadening com- 
ponent in our RRC spectral model. On the basis of the goodness 
of fit we exclude that the turbulent velocity in the recombining 
gas is greater than 3 000 km s ' . Both the temperature and emis- 
sion measure are within the error bars the same for models with 
and without a velocity broadening of 3 000 km s"' . 

If we interpret the maximum velocity broadening as due to 
Keplerian velocity at the location where the RR C is formed and 
assume a black hole mass of M = 3 x IO^Mq dMcHardv et al.l 
2004) we obtain a minimum distance for the recombining gas of 
R ~ 4x 10'^ m, or ~ 4500 Schwarzschild radii. 

From the fitted temperature of the RRC of T = 5 + 2 eV, we 
derive an ionisation parameter log^ ^ 1.6. With the ionising lu- 
minosity of the source (L = 2.8 x lO-'^ W), the minimum distance 
derived above and the definition of ^, we obtain a maximum 
density of 4x10'^ m"^. It is possible to estimate a lower limit 
to the volume of the RRC emitting source using the observed 
ionic emission mea sure of the C vi RRC (3. 9 x 10^^ m"^), a solar 
carbon abundance (lAnders & Grevess3ll98 9) and a fraction of 
60 % in the form of C vn (appropriate for log ^ = 1 .6). We find 
an emitting volume of V = 10-^^ m-*. Note that the volume might 
be larger than the emitting volume if the fil ling factor is smaller 
than unity. We should note that CLOUDY dFerland et alj|2004l) 
version 95.06 calculates the size of the clouds to be of order lO'^ 
m, thus a filling factor of less than unity is rather unlikely. We 
equate this volume to QJi^AR where AR the characteristic thick- 
ness of the emitter in the radial direction and Q. the solid angle 
sustained by the emitter. Dividing the above formula by and 
substituting the values for the volume and minimum distance we 
find: 

AR 

O — = 0.002. (1) 
R 

Taking for the distance a value of 4 x lO'"' m, instead of 4 x 
10'^ m, we derive a value of 0.02. 

Another constraint follows from the column density A^h 
through the emitting region. This can be written as A^h - nAR < 
«/?(0.002/O) (assuming the filling factor is unity, see Eq. 1). 
Taking the minimum radius and maximum density we can get 
an indication of the column density: A^h ~ 3 x 10^^(47r/Q) m"^. 

A possible origin of the RRCs is that the gas in one (or more) 
of the warm absorber components is recombining. To test this 
hypothesis, we compare the properties of the recombining gas 
with that of the warm absorber components. The large estimated 
column density for the RRC is only derived for absorption com- 
ponent 4. However for this component the ionisation parameter 
is substantially higher, making a connection between this ab- 
sorption component and the gas producing the RRCs unlikely. 
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Component 2 has a similar ionisation parameter to the RRCs but 
its column density in spectrum B and C is substantially smaller 
Thus this connection is also unlikely. Alternatively the RRCs 
could originate in the thin, ionised skin of the accretion disc. 
The width of the C vi broad emission line is similar to the max- 
imum broadening of the RRCs. Also the ionization parameters 
are similar. 

For log ^ ^ 1 .6 we expect that most of the oxygen gas will 
recombine into O vii, thus we expect a strong O vii RRC at 16.77 
A. This RRC is not observed, and we can place an upper limit 
for the emission measure of 0.6 x lO^*" m"-*. It is possible that, 
the many absorption lines due to different iron ions mask the 
presence of the O vn RRC, or that the temperature is lower than 
5eY 

lO^le et al .' (2004) detected the C vi RRC in the high flux state 
RGS spectrum; Pounds et al. (2004) detected in a different low 
state RGS spectrum the C vi and O vn RRCs. Both authors derive 
a temperature of 3 eV, consistent with our derived temperature 
of5 +2eV. 

5.5. The ionised absorber 

In our spectra, we have found evidence for the presence of four 
warm absorber components (Table [TJ. The absorption compo- 
nents are discussed in detail below. 

5.5.1. Warm absorber component 1 

The least ionised component, which has a low outflow veloc- 
ity and the smallest hydrogen column density, is not detected in 
spectrum A or B, due to the lower signal to noise in these spectra. 
The strongest Unes are due to O iv, O v and C v. 

5.5.2. Warm absorber component 2 

This component has a similarly low outflow velocity as com- 
ponent 1, but with a higher ionisation parameter, log ^=0.87. 
The strongest lines of this component are the deep ls-2p ab- 
sorption lines of Ovii, Ovi and Ov. This is the only compo- 
nent for which we measure a decrease by a factor of two in the 
ionisation parameter between spectrum C and D (see Table [T] 
FigS- This component, perhaps partly together with our com- 
po nent 1, corresponds to the low ionisation component found 
bv lKrongold et al.l (l2007l) . However, they find a column density 
that is twice as large and a somewhat higher outflow velocity of 
-500 km s'. 

5.5.3. Warm absorber component 3 

With an outflow velocity of -580 km s"', component 3 has the 
largest opacity. The dominant imprints are the Lya transitions of 
hydrogenic oxygen, nitrogen and carbon. The outflow velocity 
and ionisation stat e are similar to one of the two absorber compo- 
nents reported bv lCoUinge et al.l (1200 Ih . Note that the observed 
Doppler velocity (su m of the -h700 kms~' h eliocentric veloc- 
ity of NGC 4051, Verh eiien & Sancisil (|200l|) . and the intrinsic 
outflow velocity of ~ -600 kms"') is only +100 km s"', and 
less than 2cr different from 0. Therefore, this absorption could 
be blended with ab sorption from our Galaxy or the local group. 
[Pounds et al.l (12004) analysed the high state RGS spectrum; and 
detected an absorber with log ^-2.1 and an outflow velocity of 
~ -600 km s"' , which can probably be associated with this com- 
ponent. We associate this component with the high ionisation 



component of iKrongold et al.l ( l2007h . for which these authors 
derive a distance of 1 - 2 x 10'^ m based on time variability; but 
see Sect. 5.5.6 below. 



5.5.4. Warm absorber component 4 

The last absorber component has the highest ionisation param- 
eter, hydrogen column density and outflow velocity. The Cvi 
and Nvii Lya lines are the deepest detected. The Oviii Lya 
blends with and broadens the less blueshifted secondary Ovii 
line . The outflow velocity o f -4670 km s ' was first reported 
by van der Meer et al.] ("2003). This is one of the highest out- 
flow velocities derived from hi gh resolution spectra of Seyfert 
galaxies, that we are aware of. ICoUinee eTaP (1200 Ih observed 
for the highest ionised component in a Chandra HETG spec- 
trum of NGC 4051, an outflow velocity component of -2340 
km s ', which is not detected in our spectra. Th e high velocity 
compo nent that we detect, and the one reported bv lCoUinge et al.l 
(1200 ll) . seem to be the only ionised components detected in the 
X-rays without a UV absorber at the same outflow velocity. In 
many AGN the most highly ionised X-ray absorber should not 
produce any UV absorption lines, however usually a UV ab- 
sorber with the same outflow velocity as the X-ray absorber is 
detected. This does not seem to be the case for the -4670 km 
s"' component we detect. Both the -2340 km s"' and the -4670 
km s ' components are highly ionised an d therefore detecta ble 
with the HETG, the LETG and the RGS. lOgle et all (l200l do 
not detect these highly ionised comp onents in th e RGS spec- 
trum. Considering that neither we nor iQgle et alJ (12004) detect 
the -2340 km s"' component, and that the -4670 km s"' com- 
ponent is not present in the earlier Chandra HETG and RGS 
data, these high velocity absorbers either have a variable column 
density or are transient. 



5.5.5. Comparison between the absorber components 

The derived hydrogen column density and ionisation parameters 
for our three absorbers with an ionisation parameter log ^ <3 
are in agreement with those derived from an RGS observation, 
in a similar flux state, and given in a figure bv lOgleetalj ( l2004 . 
We confirm, using the parameters measured for spectrum C, that 
there is an increase in hydrogen column density with increasing 
ionisation parameter, as observed in this source by lOgle et akl 
(2004) and in NGC 5548 by Steenbrugge et al. (20051). 

It is also possible to constrain the opening angle of the out- 
flow. From the ionising luminosity (L = 2.8 x 10^^ W) and the 
definition of ^, we get nr^ = L/^. As we measure the outflow 
velocity v, we know that the mass outflow rate per steradian 
is m^nr^v. We list this value for each of the absorber compo- 
nents in Table 4, based on our result for spectrum C. Assuming a 
typical accretion efficiency of 10 %, the observed luminosity of 
NGC 405 1 corresponds to an accretion rate of 5 x 10""* M© yr"' . 
Assuming the mass outflow through the wind is smaller than the 
accretion onto the black hole, we can put upper limits to the solid 
angle sustained by the outflow. This assumption seems reason- 
able as the energy needed to accelerate the outflow must come 
from the infalling gas. These upper limits are small, typically 
less than 10"^ to 10"^ sr Using this same formalism, similar to 
somewhat lar ger opening angles wer e derived for the outflow in 
NGC 5548 (.Steenbrugge et al.ll2005h . 
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Table 5. Observed and calculated parameters for the four warm 
absorber components of spectrum C. 



component 


1 


2 


3 


4 


logf 


0.07 


0.87 


2.32 


3.19 


logr(K) 


4.21 


4.33 


5.14 


5.83 




2.4 X lO-*" 


3.8 X lO""* 


1.3 X lO''^ 


1.8 X 10"" 


V (kms"') 


-210 


-200 


-580 


-4670 


/Wpwr^v (Moyr"') 


1.3 


0.20 


0.021 


0.022 


a(sr) 


<3.8xlO-* 


<2.5xlO-3 


<2.4xl0-- 


<2.2xl0-' 


most abundant 










oxygen ion ; 


Ov 


Ovii 


VIII 


VIII 


T,M (s m"^) 


5.3 X 10'« 


8.8 X 10'^ 


8.5 X 10"^ 


7.4 X 10'^ 


For T, < 20 ks or 


T, > 20 ks: 








limit on r (m) 


> 9x 10'^ 


< 3 X 10'^ 


> 5 X 10"* 


> 7 X 10"* 


For T, < 570 days: 








limit on r (m) 




< 2x 10" 


< 3 X 10'* 


<4x 10'* 



5.5.6. Time variability of tine warm absorber 

Our spectral modelling of the warm absorber components is 
based upon models for photoionisation balance. Whenever the 
ionising flux of the source changes, the ionisation balance should 
adjust to the new situation, and the ionisation parameter ^ should 
change. The speed at which this happens depends on the density 
n of the absorbing material. The recombination timescale r, for 
a given ion ; is inverse proportional to t he gas density. It is given 
bv lBottorff. Korista & Shlosmanl jloOOl) to be: 

Tin = [aii— )J (2) 

where /, is the relative concentration of ion / of a given element, 
and a, is the recombination rate of ion /, and in general depends 
only on the temperature T of the gas. This formula does take 
into account the cascade into and out of ion ;. Furthermore, a 
negative value for the recombination time indicates that the ion 
is destroyed by recombination to the ion /-I. For a given pho- 
toionisation equilibrium, the temperature is a unique function of 
the ionisation parameter ^, and the ion concentrations are also 
known as a function of ^. Therefore, given a value for ^, the 
right hand side of (Eq.|2]) can be evaluated. 

The largest flux change occurs between spectrum C and D. 
The hard 2 - 10 keV power-law flux drops suddenly by a factor of 
6.1 (log p - -0.79, where we define p as the flux ratio of spec- 
trum D to spectrum C). The soft flux (taken here and later as the 
0.2 - 10 keV modified black body flux from Table 1) decreases 
by a factor of 5.4 (logp = -0.73). Thus, if the density in one 
of the components is high, we would expect log ^ to decrease by 
about 0.73, and if the density is low, we expect no change be- 
tween spectrum C and D. As the low state of spectrum D lasts 
about 20 ks and our data have insufficient S/N to resolve spec- 
trum D further in time, this gives the limit on the recombination 
timescale. 

Indeed in components 1 , 3 and 4 we see no significant change 
in log^, and definitely no change as large as -0.73 (see Table 1), 
indicating that the recombination timescale is longer than 20 ks. 
This then gives the lower limits to the distance of these compo- 
nents from the central source as listed in Table|5] In teh table the 
product nr^ is calculated from L/f . The product of recombina- 
tion timescale and density T,n is calculated from (Eq. |2]i, using 
the total recombinatio n rates of Nah ar 1999. Limits on r are 
calculated from r - •\/(nr2)(T,)/(T,«). The lower limit for com- 
ponent 3, namely 5 x lO''* m (19 light days), disagrees with the 



distance of 0.5 - 1 light days derived bv lKrongoId et al.l (l2007b 
for this component. 

iKrongold et al.l (l2007h studied the variability using the rela- 
tively low-resolution EPIC da ta. Some possib l e reas ons for the 
disagreement are given below. iKrongold et al.l (|2007|) do not al- 
low the continuum shape (i.e. photon index and temperature) to 
vary, while at least in our data the shap e of the continuum sig- 
nificantly changes (see Fig. |9l). Further, iKrongold et al.1 (I2007h 
do not include a relativistic O vm Lya line although in the total 
and high state spectra their best fit model (see their figures 2 and 
3) is consistently below the observed continuum between 17 and 
19.2 A. Similarly, the continuum is above their model between 
20 and 21.8 A, where we fit a broadened O vn line. Failure to fit 
the local continuum correctly will result in too low or too high 
column densities derived from the line depths. 

i Rrongold et alj (|2007|) included only two absorption com- 
ponents and have problem fitting some of the oxygen lines. 
Therefore, their model cannot explain the Ov absorption line 
at 22.4 A, which is easily observed in their published RGS spec- 
trum. Furthermore, the depth of the O vii line is overestimated 
in the high state, and underestimated it in the low state. Oxygen 
lines are diagnostically important because of their large oscilla- 
tor strength (i.e. strong lines) and the rather large range in ionisa- 
tion parameters they span. One can thus determine (over this ion- 
isation range) the ionisation structure without having to worry 
about abundance effects. Finally, long term differences between 
our observations and those of Krongold et al. (2007), could ex- 
plain the differences found in ionisation variability behaviour. 

iKrongold et al.l (1200 7) claim to detect variability in the ion- 
isation parameters on timescales of 2-3 ks, leading to a much 
larger derived density and thus closer location of the absorber 

In our component 2, we see a decrease of a factor of 2 
(A log ^ = -0.35 + 0.16), suggesting that the recombination time 
is indeed 20 ks or smaller This now gives an upper limit to the 
distance of 3 x 10'^ m or 110 light days (see Table |5]l. This is 
higher than, but consistent with the upper limit of 3.5 ligh t days 
for the low ionisation component of lKrongold et al.l (|2007|) . That 
the change in log ^ is not the expected value of -0.73 may have 
two reasons, which we cannot be distinguished in our present 
data. The recombination timescale might be close to our upper 
limit, so that the gas is still recombining at the end of interval D. 
Alternatively, the flux changes in the invisible EUV band might 
be smaller, and these may be more important for the low ioni- 
sation gas than the X-ray continuum (the dominant oxygen ions 
for component 2 are O vi and O v). 

Next we investigate the long term variability, by comparing 
spectrum B to spectrum C in a similar way as we did above for 
spectra D and C. The soft and hard flux of spectrum B have 
log p - -0.48 and log p - -0.33, respectively. Taking a typical 
value of -0.4, we see that all visible components in spectrum 
B (components 2-4) are in agreement with such a change in 
log^ (Alog^ = -0.38 + 0.22, -0.34 ± 0.26 and -0.27 + 0.19). 
We conclude that the recombination timescale for these compo- 
nents is likely faster than 570 days, the interval between both 
observations. This gives interesting upper limits to the distance 
of components 3 and 4 (see Table 4) of about 1 pc. We should 
note that components 2, 3 and 4 could coexist in the same lo- 
cation, but that component 1 and 2 can not coexist in the same 
location, due to the lack of variability in component 1 , even over 
570 days. 

Finally, we have an independent short term comparison by 
comparing spectrum A to spectrum B, although our discrimi- 
nating power is smaller here, because of the smaller flux jump 
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(log p — -0.13 and -0.21 for the soft and hard band, respec- 
tively). The error bars on the ionisation parameter for compo- 
nents 3 and 4 preclude a definite conclusion here; for component 
2, we see an increase in spectrum A of Alog^ - -1-0.46 + 0.32, 
which favours the no change scenario, but does not exclude a 
small response compatible with log/? = -0.13. We also note 
that, of all our four spectra, the flux during spectrum A shows 
the largest variations within the given interval (Fig. 1), so the es- 
tablishment of any (quasi)-equilibrium during that interval may 
be questioned. 

6. Summary 

In this study we have investigated the spectral properties of 
NGC 4051, observed by Chandra-LETGS on two occasions for 
a total of 180ks. Due to the flux variations within each observa- 
tion, we split each observation into 2 spectra, to be fit separately. 
In the highest signal-to-noise spectrum, C, we detect four ab- 
sorption components. Our main results are: 

- The continuum parameters do change between the different 
observations, the power-law slope becomes softer with time. 
The temperature of the modified black body increases, to first 
order, with increasing luminosity. 

- Between spectra C and D, when the flux dropped by a factor 
of 5, absorbers 1, 3 and 4 show no evidence for a change in 
ionisation parameter; however, for absorber 2 the ionisation 
parameter does decrease. 

- Our limits on the variability of the absorber suggest that 
component 2 is located within 0. 1 pc, and components 3 and 
4 are located somewhere in the range 0.02-1 pc. For compo- 
nent 1 , we only have a lower limit of 0.3 pc. These distan ces 
are different from those derived bv lKrongold et al.l (l2007h . 

- We detect one of the highest outflow velocities (- 4670 km 
s"') for an ionised absorber observed in the soft X-ray band. 
We do not confirm t he -2340 km s ' absorber detected by 
ICoUinge et al.l (1200 Ih , indicating that these high outflow ve- 
locity absorbers have variable column densities or are tran- 
sient. 

- The RRCs of C vi and C v are detected in spectrum D, af- 
ter a drop in flux by factor 5. These are unlikely due to re- 
combining gas from one of the absorber components, but are 
possibly from the accretion disc. 

- Broad emission lines of the O vii triplet and O viii Lya and 
C VI Lyor are detected. These lines show a simple Gaussian 
profile, suggesting an origin in the BLR. 

- We detect a relati vistically broadene d O vm Lya fine, similar 
to flie findings bv lOgle et all (|2004 . 
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